Strong universality and algebraic scaling in two-dimensional Ising spin glasses 
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At zero temperature, two-dimensional Ising spin glasses are known to fall into several universality 
classes. Here we consider the scaling at low but non-zero temperature and provide numerical evidence 
that r] « and v ~ 3.5 in all cases, suggesting a unique universality class. This algebraic (as opposed 
to exponential) scaling holds in particular for the ±J model, with or without dilutions and for the 
plaquette diluted model. Such a picture, associated with an exceptional behavior at T = 0, is 
consistent with a real space renormalization group approach. We also explain how the scaling of 
the specific heat is compatible with the hyperscaling prediction. 

PACS numbers: 75.10.Nr, 64.60.Fr, 75.40.Mg, 75.50.Lk 



Introduction - Critical phenomena in disordered 
systems remain a severe challenge. In fact, many prop- 
erties that are well established in pure systems (such as 
upper critical dimension, universality, etc.) are still open 
to debate when strong disorder is present. This is the 
case in particular for spin glasses (SG), i.e., magnetic 
systems incorporating both disorder and frustration Q] . 

In this work we focus on two-dimensional SG because 
(i) they are computationally tractable and (ii) their zero 
temperature (T = 0) universality classes |2( are well es- 
tablished. For our purposes there are two such T = 
classes, associated with whether or not all excitation en- 
ergies are multiple of a given quantum. The correspond- 
ing spin glass stiffness exponent 9 is zero in the first class 
(as in the ± J model) and is close to —0.28 in the second 
class (as when the couplings are Gaussian random vari- 
ables). Following standard scaling arguments 0, 0] to 
infer the finite T behavior, the thermal exponent should 
be given by v — —1/0; when 9 = 0, formally v = oo so 
that the correlation length £ should diverge faster than 
any power of inverse temperature /3 = 1/T (e.g., expo- 
nential scaling). One's expectation is then that the two 
universality classes defined at T = could have direct 
counterparts at T > 0. 

Recent numerical work 0.0 has confirmed the relation 
v = — 1/9 in the case of Gaussian couplings; similarly, all 
of the latest works 

01111013 on the ± J model ( fOT which 
9 = 0) have concluded that v = oo. We reconsider here 
this issue in greater depth, comparing models with differ- 
ent distributions of spin-spin couplings, using advanced 
Monte Carlo and partition function solvers, and applying 
sophisticated extrapolations to the infinite volume limit. 
In contrast with previous claims, we conclude that £(T) 
diverges with the same exponent v ~ 3.5 for all consid- 
ered models! Furthermore, we find that the exponent rj 
also seems to be the same regardless of the details of the 
couplings. We are thus in a situation where the T > 



properties fall into just one universality class, whereas 
two classes arise when considering the T = properties 
(where the quantized nature of energies can be relevant). 
A simple interpretation of this unusual phenomenon is 
that at T = one has an "additional" fixed point (re- 
alized only when energies are quantized) but that this 
fixed point is irrelevant for the critical properties defined 
just above T c . We have also investigated the specific heat 
singularities of our models. Here we find a behavior that 
is asymptotically compatible with T c = hyperscaling. 

Models, computational methods and, analysis - We 
consider the 2D Edwards- Anderson model defined by the 
Hamiltonian 
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where the sum is over all the nearest neighbor pairs of 
a 2D square lattice of length L with periodic bound- 
ary conditions. The a x = ±1 are Ising spins. We take 
the couplings J xy to be independent random variables, 
chosen from several distributions in order to check their 
effect on the critical behavior. In particular, we shall 
consider: (1) the ± J model where J xy = ±1 with equal 
probability; (2) the diluted ± J model obtained from the 
undiluted case by setting a fraction (1 — p) of the cou- 
plings to zero; (3) the "irrational model" where J xy = ±1 
or ±G with equal probability for these 4 values, G being 
the golden mean, G = (1 - \/5)/2 ~ -0.618; (4) A distri- 
bution with i of the couplings equal to ± 1 and one half 
equal to ±|: we call it "Gap 1/4". 

We tackled these systems using two complementary ap- 
proaches. First, we applied cluster replica Monte Carlo 
along with parallel tempering and numerous algorithmic 
optimizations H 0, ^| ; such an approach can provide 
high quality estimates of configurational averages even 
for large lattices. Second, we used exact partition func- 
tion computations 0] which give us the value of the free 
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energy and specific heat for individual samples at arbi- 
trary temperatures. In all cases, for each lattice size L, 
we used a large number of disorder realizations for each 
of the distributions mentioned previously; typically, we 
used thousands of such samples. In spite of the large 
lattices sizes we are able to handle, it is appropriate to 
use sophisticated techniques for analyzing the effects of 
finite L. Starting with a suitably defined finite-volume 
correlation length £(T,L) (see later) and a long-range 
observable 0(T,L) such as the spin glass susceptibility 
Xsg(T,L), finite size scaling (FSS) theory predicts 



Q(T, sL) 
0(T, L) 



(t{T,L)/L; 
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(2) 



Here, Fq is the FSS function and s > 1 is a scale factor. 
Eq. J2J) is an excellent starting point for investigations of 
the FSS behavior, as it involves only finite- volume quan- 
tities taken from a pair of systems with sizes L and sL 
at a given T. The knowledge of the scaling functions 
Fq (where O is in our case xsg) an d F$ (where O is £) 
allows us to extract information on the critical behav- 
ior using an infinite volume extrapolation 

This 

technique works with data strictly above T c and hence is 
well suited to our case for which T c = 0. In essence it 
uses and Fq to obtain the thermodynamic limit of O 
using an iterative procedure in which the pair £ and O is 
scaled up from L — > sL — > s 2 L — > . . . — > oo as described 
in ^]|. Detailed knowledge of the critical behavior can 
then be obtained from appropriate fits to the extrapo- 
lated data. To show that this relatively sophisticated 
technique works well in our system, we display in Fig. ^ 
the FSS function F$ determined from these procedures, 
for several of our models. The excellent data collapse 
at small values of £/L (£,/L < 0.45) (already noticed in 
[r?l llijl validates the FSS framework. Furthermore we 
see that these FSS functions are independent of the dis- 
tribution of couplings J xy ; that is precisely what should 
transpire if there is a single universality class in our sys- 
tem. The region of the excellent data collapse increases 
with the size of the system suggesting strongly a sin- 
gle limiting FSS function. The diluted and the Gap 1/4 
model are the most important pieces of evidence on which 
our conclusions rest. The diluted model has smaller fi- 
nite size effects than the undiluted one: this allows us to 
go to very reliable extrapolations. We have been able to 
thermalize systems up to L — 64 (with 2000 samples). 
The spin-glass susceptibility and the exponent i] - 
In our Monte Carlo simulations we follow two repli- 
cas from which we can measure the local spin overlaps 
Qxy = cTxyCTxy (the superscript is the replica index). Of 
interest is the two site correlation function of these over- 
laps; we define in fact the wave- vector-dependent spin- 
glass susceptibility via 



,ik'(r2—ri) 



r 2 /Tjav 



(3) 
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FIG. 1: The FSS function F 6 . 



is then defined 
= 0). (We denote 
)t and the 



The usual spin-glass susceptibility 
through X sg(T,L) = L 2 (q 2 ) = X sc{k 
the thermal average at temperature T by 
average over the disorder realizations by [. . .] av .) 

The spin-glass susceptibility is a measure of the corre- 
lation volume of our system; using the standard form for 
a correlation function near a critical point, 



G(r) 



r 
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we see that in D = 2 as T — > XSG should behave as 
XSG ~ £ 2 ~ ,? . To check this, we use our infinite volume ex- 
trapolations for both £(T) and xsg(T), and display one 
versus the other in Fig.|2 It is remarkable that the data 
for our models fall on the same curve, indicating in par- 
ticular that they all have the same exponent 77, a highly 
surprising fact if there were two universality classes. Fits 
of this curve lead to values of 77 that are very small, be- 
tween and 0.1, strongly suggestive of 77 = 0; note that 
one expects the models with continuous couplings to have 
i] = exactly. 

The correlation length and the exponent v — The sec- 
ond moment of the overlap correlation function can be 
identified with the square of the correlation length. In 
practice, we define this length via 
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where xsg is defined in Eq. ©. Note that k m i n = 2n/L 
is the modulus of the smallest non-zero wave-vector al- 
lowed by periodic boundary conditions. We use our anal- 
ysis methods to extract the large L limit of £(T, L). Hav- 
ing done so for our different models, we find that their 
temperature dependence is not so different. This is illus- 
trated in Fig. 13 we also see that the different models lead 
to rather parallel asymptotes, suggesting that the expo- 
nent v is the same for all of them. Analogous plots on 
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a semi-log scale do not give evidence for an exponential 
scaling of £, even for the ± J model. Finally, performing 
power fits, we find values of v close to 3.5. Again, this 
goes in the direction of a single universality class. 



1 Undiluted +/-J 

. Diluted +/-J, p = 0.875 
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♦ Irrational 

-- T| =0 




FIG. 2: Parameter-free plot of Xsg(T, oo) versus £(T, oo) for 
different coupling distributions. The line is for 77 = 0. 

The specific heat and the exponent a - We mea- 
sured the specific heat density c v both via Monte Carlo 
and via the exact partition function techniques. In con- 
trast to what was observed for x an d £, the ±J model 
and its diluted versions differ from the continuous mod- 
els when one considers the very low T behavior of c v . 
In the first class of models, c v decreases fast as one ap- 
proaches T — 0, while in the second class, c v looks linear 
in T. According to hyperscaling, when T c = 0, one has 
the relation 191 



c v {T) ~ T~ 



with 



o 



-2v 



(6) 



We saw that v w 3.5, so we expect a w —7. The decep- 
tive linear behavior in the continuous model can be easily 
understood. In the ground state of a continuous model, 
the local field on a spin is a continuous variable with a fi- 
nite density at 0. Thus at any low temperature, a fraction 
proportional to T of the system's spins can be thermally 
activated, thereby leading to a linear dependence of c v 
in T. The point is that the exponent a describes the 
non-analytic scaling of c Vl but one can also have regular 
parts, namely linear, quadratic etc. in T . Since —a is 
large, the non-analytic part is sub-dominant compared 
to the analytic contributions and thus numerically invis- 
ible. The conclusion is that the continuous models make 
it impossible in practice to estimate a. Fortunately, the 
situation is completely different for the ± J model and its 
variants. There the leading analytic contribution to c v 
is of the form T~ 2 exp(— 4 J/T) because the local field on 
any spin is a multiple of 2 J. This is far smaller than the 
predicted non-analytic term (a power of T), so we are 
able to use these models to estimate a. In Fig. 0]we 
show the specific heat of the ± J model versus T. In the 
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FIG. 3: £(/?, oo) as a function of j3. 
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FIG. 4: Specific heat of the ±J model versus T in a log-log 
plot. Inset: plot of -Tln(T 2 c„) versus T. 



inset we show the semi-logarithmic plot of c v to test for 
exponential scaling. The data are not compatible with an 
exponential scaling (the L = 80 data are crucial to reach 
this conclusion), while an asymptotic power law behavior 
is consistent with the data and with the hyperscaling of 
Eq. ©. 2D disordered systems with different choices of 
a (frustrated) disorder behave in the same way 20], mak- 
ing our claim about a strong universality even stronger. 

A renormalization group justification - In the pic- 
ture we are proposing the thermal properties (T > 0) of 
2D SG fall into just one universality class, in contrast to 
the properties defined at exactly T = 0. Furthermore, 
the relation v = —1/0 holds if we take 6 from the T = 
class of "continuous models" (where 8 is non-zero). We 
are thus faced with a situation where continuous models 
behave as expected, but the ± J models have (i) a sep- 
arate (exceptional) behavior at T = 0; (ii) fall into the 
continuous class when T > 0. To justify theoretically 
this scenario, we appeal to a real space renormalization 
group (RG) framework. Consider the Migdal-Kadanoff 
(MK) construction for spin glasses [2]]]. These hierar- 
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chical lattices can be generated by recursively applying 
an "inflation" operation whereby each edge is expanded 
into a cell having b parallel branches of length s. Renor- 
malization on these lattices can be done exactly. For 
the purpose of this discussion, consider the case where 
b = s = 3. When the distribution of couplings is con- 
tinuous, one has 8 « —0.278, while when using the ±J 
distribution one has 8 — 0. Indeed there are two fixed 
point distributions of the couplings: a generic (continu- 
ous) one for which the energy scale goes to zero (8 < 0) 
and another one which is very special, being a sum of 
delta functions on odd integers (and leading to 9 = 0). 

Now what happens when we turn on the temperature 
in this real space RG framework? At very low T one is 
very close to the critical manifold so the renormalization 
will first flow toward one of the T = fixed points; the 
corresponding fixed point determines the critical expo- 
nents of the system. In the case of the ± J model, the 
initial distribution starts by being nearly concentrated 
on odd integers. However, as we increase the lattice size, 
this initial distribution will be less and less of that form: 
the iterations renormalize the free energy, F = E — TS, 
and since S fluctuates, we loose the quantization property 
rapidly with size. The flow is thus toward the "continu- 
ous" distribution fixed point for which 8 < 0. In Fig. 
we display the specific heat as a function of T for this ±J 
model; as expected, the behavior indicates a power scal- 
ing of c v totally compatible with the value of a obtained 
from the continuous distribution fixed point. 
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FIG. 5: Specific heat for the ± J model on MK lattices versus 
T in a log-log display. G is the number of recursions of the 
MK construction. Inset: same data but for a semi-log display, 
which should be straight in the case of exponential scaling. 

Summary and conclusions — In this study of two di- 
mensional SG, we examined the effect of the underlying 
distribution of couplings on the critical thermodynamics. 
We found several remarkable behaviors: (1) n is model 
insensitive and very close to zero in all models. (2) the 
divergence of £ does not seem to depend on the model, 
and the corresponding exponent is given by v = —1/8 



as long as one uses the 8 of the continuous models. (3) 
The specific heat exponent a can be measured using the 
discrete models, whereas in the continuous models this 
singularity is sub-dominant and thus numerically inac- 
cessible. Hyperscaling is consistent with the data. 

From our observations we conclude that thermal prop- 
erties of 2D SG fall into a single universality class. This is 
in contrast to what happens exactly at T = 0: there, an 
additional fixed point of the RG arises, allowing different 
scaling in the continuous and in the ± J type models. 
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